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Abstract 

We explore several models of QCD evolution equations simplified by considering only the 
rapidity dependence of dipole scattering amplitudes, while provisionally neglecting their de¬ 
pendence on transverse coordinates. Our main focus is on the equations that include the 
processes of pomeron splittings. We examine the algebraic structures of the governing equa¬ 
tion hierarchies, as well as the asymptotic behavior of their solutions in the large-rapidity 
limit. 


1 Permanent address. 



1 Introduction 


Parton saturation and the unitarity bound are among the interesting, yet still not fully resolved 
problems of the high energy QCD. There are two equivalent descriptions of parton saturation 
effects: the Balitsky hierarchy, . 1 ] and the JIMWLK functional equation 0 - The Balitsky hierarchy 
consists of an infinite set of coupled evolution equations for the operators involving Wilson lines. 
They describe quark and gluon scattering off a target, and generally are applicable in the high- 
energy regime. On the other hand, the JIMWLK equation is a functional equation which describes 
the evolution of the target in the limit of small x. When the JIMWLK hamiltonian acts onto a 
dipole operator consisting of Wilson lines in the fundamental representation, the Balitsky hierarchy 
is reproduced. In the limit of negligible correlations within the target, the first equation of the 
Balitsky hierarchy decouples, and becomes equivalent to the nonlinear Kovchegov equation 0. 
The latter has been derived within the dipole formalism 0 by taking into the account multiple 
scatterings between the dipoles in the projectile wavefunction and the target. 

A common feature of the above approaches is that they involve triple interactions of BFKL 
pomerons, [5] and thus describe the evolution of the system as a nonlinear process. However, it has 
been recently recognized that the Balitsky-JIMWLK equations incorporate only pomeron mergings, 
but not pomeron splittings 0. In consequence, in the course of the evolution (understood as a 
process in the rapidity parameter) the number of pomerons can only decrease. Therefore the above 
equations ignore contributions from fluctuations from any intermediate pomeron loops. However, 
it is well known that such terms are relevant and can be substantial, therefore they should not be 
neglected 00. Accordingly, considerable efforts have been lately undertaken towards a better 
understanding of the role of the fluctuations, Hum by including the missing pomeron splittings 
in the Balitsky-JIMWLK equations mmmm 

Let us first consider the evolution equations as proposed in 0. They describe the evolution of 
(k) 

the amplitudes Ty ' for the scattering of k quark-antiquark dipoles off a target in the large number 
of colors N c —> oo regime. The rapidity Y = In 1/a; substitutes for the (imaginary) time parameter 
of the evolution. The first two equations of the infinite hierarchy read 


M(x,y,z) ® T| 1} (x,y) - Ki(x, y, z)t£ 2) (x, z; z, y) , 

M(x i,yi,z) ® Ty- 2) (xi,yi;X 2 ,y 2 ) + (1) 

+ K 2 (xi,yi,x 2 ,y 2 ,z) ® T| 1} (xi,yi) - Ki(xi,yi, z):Z\ ( - 3) (...) . 

(k) 

The arguments of the dipole amplitudes Ty ; are two-dimensional vectors x and y that give the 
location of the dipoles in the position space. The value of the strong coupling constant a s = as ^ c 
remains fixed, consistent with the leading lnl/a: approximation. The functions and K 2 

describe the BFKL kernel and the Pomeron vertices in the dipole picture, respectively. Their 
exact forms can be found in 0. Note that the first equation is coupled through Ty' 1 to the second 
one, the second one is coupled through Ty ; to the third one, and so on. This way the infinite 
hierarchy is formed. 

Describing the exact behavior of the solutions to this hierarchy poses a nontrivial task. Even in 
the case of the closed nonlinear Kovchegov equation, the exact analytical solution is not available, 


ai Y {x,y) 

dY 


dTy ] (xi, yi; x 2 , y 2 ) 


dY 


d 2 2 
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although one can gain some partial insight about its behavior from analytic properties of various 
approximate cases mmm, or numerical simulations ( see for example El). 

In this contribution we study the simplified model of the infinite hierarchy 0, in which we 
ignore the dependence on transverse coordinates and retain only the dependence on the evolution 
variable Y. This is of course a very restrictive approximation, yet nevertheless it will enable us 
to explore the structure of the hierarchies, the broad properties of the solutions, as well as their 
rapidity dependence. In particular, we will demonstrate how to construct in a relatively simple 
way the general solutions to the infinite hierarchy with, and without fluctuations. We will show 
that in one particular model with fluctuations included, the solutions can be written in a compact 
form with the differential operators expressible as Fibonacci polynomials. We will also show that 
in a model that is close to the new hierarchy m the fluctuations are significant in the regime of 
small amplitudes. We point out however, that in this case the physically interesting solutions do 
not possess a common limiting point and that they are valid only as long as k < 1/oq, where k 
denotes the number of dipoles. 


2 Model I without fluctuations 


Let us consider a simplified model for the infinite hierarchy of equations 


dx k jt) 
dt 


x k (t) - x k +i(t) , 


( 2 ) 


where x k {t ) is the amplitude for the scattering of k dipoles off the target. Here it depends only on 
the evolution variable t. The first term on the right-hand side is responsible for the free propagation, 
while the second term describes the process responsible for pomeron merging. 

Let us rewrite the equation © in a matrix form, 


j t Mt) = M X(i), 


( 3 ) 


where 


X(t) 


/ Xi(t) \ 
X2 (t) 

X3 (t) 


X n (t ) 


v • 


is the vector of the solutions for 
form 


M = 


the amplitudes 

/ 1 — 1 0 

0 1 -1 

0 0 1 

0 0 0 

0 0 0 

0 0 0 


X k(t). 

0 

0 

-1 

1 

0 

0 


In the present model the matrix M has the 

0 0 ... \ 

0 0 ... 

0 0 ... 

-1 0 ... 

1 -1 ... 

0 1 ... 


V 


/ 
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Note that every finite matrix M* n) that corresponds to the truncated version of m consisting 
of just n equations has the determinant detM^ = 1 and has n degenerate eigenvalues Ai = • • • = 
A n = 1. However, in order to solve the complete set ©, it is required to find the eigenvectors of 
the full infinite matrix M. It is straightforward to notice that the vector 

\ 




with A = 1. This defines a particular solution of the infinite hierarchy © 

xi (f) = e t , 

Xk(t) = 0 , k > 1 , 

which corresponds to the free propagation of a single particle analogous to the BFKL pomeron. 
To find a more general set of solutions we postulate that 

xi (t) = e 4 /(f) , ( 6 ) 

where /(f) represents an arbitrary function. By recursive use of the equations (J2J we find that 

x k+1 (t) = (-l)V^, fc> 0 , (7) 

jO £ 

where = /■ So far, the only restriction on the function /(t) is that it belongs to C°° . These 
solutions can be given a more compact form with help of a generating functional Q for x k 

OO 

Q[u] = ^2 uk Xk +1 , ( 8 ) 

k—0 

or exactly, 

S[u] = e 4 —(9) 

1 + U dt 

where it is understood that the differential operator acts only on the function /(f) which in fact 
determines the specific solution. For example, the choice of /(f) = e -4 corresponds to constant 
solutions 

Xl = X2 = • ' ' = Xk = ' ' ' = 1 • (10) 

In order to find a nontrivial solution, let us provisionally assume that 

x 2 = 'yxiXh ( 11 ) 

where 7 is a non-zero constant. In the case of 7 ^ 1 , it can be interpreted as a source of non-zero 
correlations DU EE] From the assumption GH> it follows that the function /(f) satisfies the 
differential equation 

I = ‘ ’ < 12 > 


(4) 

(5) 


( 1 
0 
0 


V = 


solves the eigenvalue problem 


M V = A'' 
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whose solution is 


1 1 


(13) 


fit) = 

Using this solution in 0 we readily find that 


Xi(t) = 

1 

e 4 

7 f 

4 + C ’ 

X2 (t) = 

1 

( e 4 \ 

7 

Kc + c) 

x 3 (t) = 

1 

2 e 34 

7 

_(e 4 + C) 

X 4 (t) = 

1 

CO 

1 _ 

7 

(e 4 + C) 


7 e 4 + C 


(e 4 + Cf 
6 e 34 

(e 4 + Cf 


(e 4 + C) 2 


(14) 


Interestingly, in the t —> oo limit, all the amplitudes x k —■> 1 / 7 . In order to satisfy the unitarity 
bound for the amplitudes, it is required that Xk < 1 and, accordingly, 7 > 1 . 


3 Model II for the hierarchy with fluctuations 


Let us now modify the hierarchy by enriching it with terms responsible for particle splitting [SOU 

Ej 


(with xq = 0). Again, 

where now 


or symbolically, 


dxkjt) 

dt 


M = 


= X k (t) - X k +i(t) + X k -i{t), k> 1 , 


|x(i) = MX(f) , 

/ 1 -1 0 0 0 0 ... \ 

11-1 0 0 0 ... 

0 1 1-10 0 ... 

0 0 1 1 -1 0 ... 

0 0 0 1 1 -1 ... 

0 0 0 0 1 1 ... 

V / 

M = 1 - A + A t , 


(15) 

(16) 


(17) 


where Aij = <5, y —1 - The elements of A are responsible for transitions from k to k — 1 particles, 
while the elements of A T generate transitions from k — 1 to k particles. It is straightforward to 
construct a generalization to a more complicated structure with k-to-n particle transitions with 
arbitrary k 7 n once the corresponding vertices are known. 

The problem now becomes somewhat less degenerate since the finite (n x n) submatrices M^ n) 
of the truncated hierarchy have different eigenvalues. We also remark that in the present case the 
determinant det M is infinite. 
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The fact that the matrix M has only a finite number of non-zero terms in any given row or 
column makes it easy to construct an eigenvector. For example, the eigenvector to the eigenvalue 
A = 1 of the problem 

MW = AW , (18) 

1 
0 
1 
0 
1 
0 

V • 

We note that this vector is not normalizable; wtw = oo. Nevertheless, we shall adopt it as a 
particular solution of CU, whose exact form is 

X 2 i-i{t) = e*, 

*2i (*) = 0, 1>1. (19) 

Again, one can construct a generalized solution by taking an arbitrary function f(t), postulating 
that x\ (t) = e t f{t), and using the hierarchy equations to generate the remaining Xk(t). As before, 
the specific form of f{t) is dictated by the type of the solutions sought. In the general case, the 
first six amplitudes are 


is 


W = 


*i (t) = e‘/ , 

X 2 (t) = -e‘/ (1) , 

*3(t) = e‘(/ + / (2) ) , 

* 4 (t) = -e*(2/ (1) +/ (3) ) , 

*b(*) = e* (/ + 3/W + /W) , 

*eW = -e 4 (3/W+4/( 3 ) + /( 5 )) , 


( 20 ) 


where 

The perusal of the coefficients of the derivatives reveals that they are the same as in the 
Fibonacci polynomials, therefore we can rewrite the solution in a more compact form, 

x k+l (t) = (-1 ) fc e 4 ^[A]/(i), fc > 0 , (21) 

where Tk[u] are the Fibonacci polynomials (for the case of Tq[u] = 1). The generating functional 
does not become much more complicated, and it reads 

G[u) = e t — l -- f(t ) . (22) 

1 + u Tt ~ u 
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If we again assume the approximate factorization EJ then the corresponding particular solution 


Xi(t) 

27 (i) 
X3 (t) 

Xi(t) 

x 5 {t) 


C 


1 e* 

7 e* + C ’ 

1 / ' 2 
7 
1 
7 
1 
7 
1 
7 


24 


e* + (7 
e* 

e‘ + C 
e* 

e* + C 


e* + C 


e* + C 


- 6 


e* + C 


+ 3 


e* + C 


- 36 


e* + C 


20 


e* + <7 


-4 


( 23 ) 


C 


+ 


c 


Unlike in the previous hierarchy, the solutions Xk do not approach the common l/y limit at t —> 00. 
Instead, (for simplicity we take 7 = 1) the solution is 


27(00) 
2:3(00) 
2:4(00) 
2:5(00) 
Xk( 00) 


27(00) = 1 , 

2:2(00) + 27(00) = 2 , 
2:3(00) + 2:2(00) = 3 , 
2:4(00) + 2:3(00) = 5 , 
27-1(00) +x k - 2 {oo). 


( 24 ) 


which we recognize to be the Fibonacci series. One of the properties of the Fibonacci series is that 
it exhibits the exponential behavior, 


F n 



( 25 ) 


where a — 1.6 and [... ]i nt stands for the nearest integer value. In short, in this simple model, 
when the amplitudes 27 and 27 reach the limit 


27,27 = 1 , 

then the dipole amplitudes of higher-order in k must exceed this limit. However, in order to be 
physically meaningful, the spectrum of solutions should be bounded, 

0 < 27 < 1, \/k > 1 , (26) 

which generally imposes an implicit condition on the behavior of f(t) in ( 11211 ) . 

Unfortunately, this inconsistency is fatal. The condition of boundedness of 27 requires that the 
derivatives vanish as t — > 00 . In this limit the hierarchy simplifies to 

27+1 — 27 -f- 27_i, (^7) 

which, irrespectively of the form of the function /(£), is precisely the diverging Fibonacci series 

(T?51) . 
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4 Model III for the Balitsky hierarchy 


In this section we consider the model for the simplified Balitsky hierarchy [T] 

dx k {t ) 


dt 


= k [x k (t) - x k +i(t)] 


(28) 


(compare with the equation In this case the matrix M is still triangular, although the 

corresponding determinant is infinite. For the truncated hierarchy it behaves as detM* n) = nl. 
The eigenvector corresponding to the eigenvalue A = 1 is the same as in the former case, 

f 1 \ 

0 
0 

V = 


V ■ 

Therefore we can reaply the procedure developed in the last sections to obtain 

A 


Xi(t) 

x 2 (t) 


(29) 


x 3 (t) 
X4 (t) 


= e - 


(; 5 (/“>+ /<’>)-/“>) , 

f m + \ ( f m + / <2) ) + j (/<■> + / <2 > + I (-/ (I) - 2/<=> 


The factorization ansatz x 2 = 7 x 1 x 1 leads to 


Xk(t) = ~ 
7 


C 


(30) 


which all converge to the the same limiting point Xfc(oo) = 1 when 7 = 1 . We recognize that the 
above are indeed the mean-field solutions of the model found in mm- The first six solutions of 
this sequence are illustrated in Fig. Q ■ 


5 Model IV for the Balitsky hierarchy with fluctuations 


Let us proceed to the next model in which we add also the fluctuations terms 2 to iPl) 
dx k (t) 


dt 


= k [x k (t) - Xfc+i(t)] +a 2 s k(k- 1) Xk_i(t) . 


(31) 


The form of the coefficient a 2 s k{k — 1) is motivated by the structure of ©, see JJ. for details. In 
this case the matrix elements are 


M-mn TYl8mn OL s 7fl{7Tl 1 )—In • 


(32) 


2 The solutions to the hierarchy with Pomeron loops have been recently investigated in im 









Figure 1: The factorizable solutions to the toy model for the Balitsky hierarchy iPl). From left to 
right: x\,..., x$. The constant of integration has been set at C = 10, and the coefficient 7 = 1 . 


The recurrence relation for the static solution 


dx 


(s) 


dt 


= 0 is 


°k +1 


= 4 S) + a K k - > 


which for large values of k leads to 


(33) 


4^ ~ ( a s Vk) k . 


(34) 


This is even stronger divergence than in the case of the Fibonacci series. As long as k is not 
too large, k <C l/cc 2 , the magnitudes of these static solutions remain close to each other, but for 
k's larger than that, they quickly diverge. 

One can solve the recurrence relation by means of generating functional which in this case 
we define as 

sm = E»"4. < 35 > 

n= 1 

(compare the definition JHJl). The differential equation for the functional reads 


d 2 g 

du 2 


dG_ 

du 


(1 + a 2 s u) 


and has the solution 

g\u] = f dye v+olsV ! 2 , 

J 0 


which can be expressed in terms of a confluent hypergeometric function \F\ 


g[u\ 



1F1 


(1 3 (l + a^) 2 \ 

V2’2’ 2^ J 


1 F 1 


1 3 1 V 


(36) 


The physically interesting solutions need to be positive-definite functions. We find that the solution 
X\ (t) from (THll gives positive definite solutions Xk(t) also in the case of the hierarchy (131 [l . The 
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Figure 2: Left plot: the first six solutions to the toy model hierarchy (ZD- The curves from the 
left to the right: x\, ..., x$. The coupling has been set to a s = 0.2, and the coefficient 7 = 1. The 

constant of integration wass set at C = 10. Right: The ratio of the solutions 01 to the factorized 

( f) 

solutions of the Balitsky hierarchy; R k = Xk (t)/x k '(t). The constant line at R k = 1 is the ratio for 
x\/x{ and Xz/x^- The four curves from bottom to top are xs/x^\ xa/x^\ x^/x^ and xq/x^\ 


first four solutions are 


Xl (t) 
x-i (t) 
X3 (*) 
XA (t) 


1 e 4 


7 e 4 + C 
1 / ^ ' 2 


7 \ e* + C, 

1 / ' 3 


7 \e t + C; 

1 / ' 4 


7 \ e* + C 


The general form of the solutions is 

x 
(f) 

where x\ (t) are factorizable solutions 


(37) 


7 e‘ + C ’ 
al e t (9e t + 2 C) 
37 (e 4 + C) 2 


Xk{t ) = x[ f \t) + Ax^\t) 


(38) 


x[ f \t) = (xi(t)) k , 

given by OT . The corrections Axjj l) ( t ) are nonfactorizable and suppressed by a factor of at least 
a 2 s . They do not vanish at t —> 00 , and contribute non-negligible corrections to the magnitudes of 
the limiting points. We illustrate these particular solutions (El on the right plot in Fig. [2j where 
the differences between the positions of the limiting points Xk( 00 ) are apparent. It is interesting 
to look in better detail at the magnitudes of the new nonfactorizable terms. For this purpose, in 

the right plot in Fig. [2 we show the ratio of the solutions Xk/x): as a function of time t for the 

( f) ( f) 

first six solutions. The solid line at 1 corresponds to the constant ratios of xi/x\ J and X 2 /X 2 ■ 
Beginning from 23 , the non-factorizable terms are nonzero, and they are most important at very 
low values of t, where the amplitudes Xk are small, as expected ®ma. We also consider an 
alternative model with an additional negative term proportional to Xk , 


dx k (t) 

dt 


k [a : k {t) - x k +i (t)] +a 2 s k(k- 1 ) [x k -\ (t) - x k (t)] . 


(39) 
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Figure 3: Left: the first six solutions to the toy model hierarchy cu. The curves from the left 
to the right are The coupling has been set to a s = 0.2 and the coefficient 7 = 1. 

The constant of integration set at C = 10. Right: The ratio of the solutions El to the factorized 

( f) 

solutions of the Balitsky hierarchy; Rj- = Xk(t)/x k '(t). The constant line at Rk = 1 is the ratio for 
x\/x{ and x 2 /x 2 - The four curves from bottom to top are x 3 /x^\xi/x^\xs/x^ and xq/x^\ 

This particular form of the hierarchy is inspired by m- The role of the additional term is to secure 
the conservation of probability in the pomeron splitting processes. We must however remark that 
this term is formally proportional to the BFKL kernel generated by the not yet available next- 
to-next-to-leading order logarithmic approximation. With the negative term included, the static 
solution to El reads 

44 = 4 S) + ( fc -!) (4-i - 4 s) ) > ( 4 °) 

and has a nice property that 

= • • • = x^ = • • • = const , (41) 

i.e. it predicts the same limiting point for all particular solutions xjt, just like in the model El- 
Now the particular solutions that satisfy the condition X 2 = 7 x 1 X 1 have the form 

. . 1 e* 


(42) 

a 2 Ce t 
7 ( e t + C) 2 

a 2 2e t C(C 2 + 6Ce 4 + 5e 24 ) a\ 26V 

— 3(e 4 + C) 4 T (e 4 + C) 2 ' 

As before, we can express the solutions in the general form Xfc(f) = x^\t) + A x^\t). In this 
case the nonfactorizable terms vanish at least as e -4 for t —> 00, which guarantees the convergence 
to the same static solution GD- This behavior is illustrated in Fig. El The left plot in Fig. El 
shows the first six solutions El- The right plot shows the ratio of the solutions Xk/x\f\ which 
behaves similarily as in the case of the hierarchy ED- The aditional terms ~ a 2 responsible for 
the pomeron loops are most significant in the low t regime, where it<l. 


xi(i) 

x 2 (i) 

X3 (t) 
X4 (t) 


7 e 4 + C 
If e 4 


7 \e 4 + C/ ’ 


If e 


t \ 3 


7 \e 4 + C / 
If e 4 x 4 


7 Ve 4 + 6' 
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Conclusions 


We have shown that the general features of the asymptotic solutions of the hierarchies for the 
evolution equations in QCD can be deduced from simplified models that keep only the rapidity 
dependence while neglect the dependence on the transverse degrees of freedom. We presented the 
general method for finding the solutions to these types of hierarchies. It is generally expected that 
on the grounds of the unitarity or ’black disc’ conditions the solutions ought to tend to a common 
limit. However we have found that in some cases of hierarchies with fluctuations the physically 
relevant solutions may not tend to one universal asymptotic solution. In some models they may 
even form highly divergent series. We have shown that the fluctuations have the largest impact in 
the region of small values of amplitudes, irrespectively of the model considered. 
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